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Prospects for Lunar Satellite Detection of Radio Pulses from 
Ultrahigh Energy Neutrinos Interacting with the Moon 

O. Stal, 1 - 2 '0 J- E. S. Bergman, 1 B. Thide, 1 - 3 L. K. S. Daldorff, 1 - 3 and G. Ingelman 2 
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The Moon provides a huge effective detector volume for ultrahigh energy cosmic neutrinos, which generate 
coherent radio pulses in the lunar surface layer due to the Askaryan effect. In light of presently considered 
lunar missions, we propose radio measurements from a Moon-orbiting satellite. First systematic Monte Carlo 
simulations demonstrate the detectability of Askaryan pulses from neutrinos with energies above 10 eV, i.e. 
near and above the interesting GZK limit, at the very low fluxes predicted in different scenarios. 
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Ultrahigh energy cosmic rays (UHECR) provide impor- 
tant information on extreme astrophysical and cosmological 
processes. Above the Greisen-Zatsepin-Kuzmin (GZK) limit 
~5x 10 19 eV the universe should be opaque over intergalactic 
scales to protons due to py — > N7t where y is a 2.7 K cosmic 
microwave background photon [1]. The subsequent decays 
7T + — > [J. + "Vu and (i — > WeVfi gi ve an isotropic flux of "GZK 
neutrinos" ||2J,|3f] which, if detected, should resolve |4J]the ap- 
parent contradiction implied by recent observations |5J, |g] of 
cosmic rays with energies above the GZK limit. 

Weakly interacting ultrahigh energy neutrinos (UHECv) 
propagate unaffected over cosmic distances, so their arrival 
directions point back to the original sources. Detection of 
neutrinos with energies well beyond the GZK limit has also 
been suggested as a method to test cosmology through the 
Z-burst process [7], in which the highest energy cosmic rays 
would be produced following resonant interactions vv — ► Z° 
of UHECv with cosmological relic neutrinos. In addition to 
these mechanisms, UHECR and UHECv could be produced 
in the decay of super-heavy relic "X" particles originating e.g. 
from topological defects (TD) associated with Grand Unified 
Theories Q8Q. Recent neutrino flux limits, in particular that 
from ANITA-lite |9fl, constrain these models and allow only a 
percent-level contribution to the UHECR from Z-bursts. 

No UHECv has been observed with existing neutrino tele- 
scopes. The next generation km 3 -sized optical detector Ice- 
Cube II 1 Of! is optimized for 10 12 — 10 17 eV neutrinos. At 
even higher energies, it seems more promising to explore 
the Askaryan effect 111 111 of Vavilov-Cerenkov (VC) radio 
pulse emission from a charged particle shower in a dielectric 
medium. The shower, initiated by the interaction of a high en- 
ergy particle, produces via secondary scattering in the medium 
a net charge excess which radiates coherently for wavelengths 
longer than the shower dimension. Hence, the power radiated 
in radio scales quadratically with the initial particle energy 
and not linearly as in the optical region. 

The extremely low flux of UHEC"v necessitates a huge 
effective target. This can be achieved by the detection of 
Askaryan pulses from neutrinos interacting with the upper 
layer (regolith) of the Moon, which has been considered for 



Earth-based radio measurements 1 12lll3lll4ll5lllo1.il 711 . and 
is currently proposed for a lunar satellite mission 11811 . Here, 
we report Monte Carlo (MC) simulation results on Askaryan 
pulse detectability with instruments on a satellite orbiting the 
Moon. 

The primary neutrino-nucleon cross sections for deep in- 
elastic charged current (CC) and neutral current (NC) inter- 
actions are well-known and can be extrapolated to higher en- 
ergies 11911 . Both CC and NC interactions initiate hadronic 
showers carrying E s = yE v of the total neutrino energy E v , 
with (y) = 0.2 at the very highest energies 11911 . As discussed 
below, we only consider hadronic showers at this stage. 

The Askaryan radio emissions from showers in different di- 
electric media have been parametrized 120il2U 12211 and the re- 
sults for Silica — the main constituent of the lunar regolith — 



have been validated in accelerator experiments 123l 12411 . In 
the form given by 11711 . the spectral flux density F in Jansky 
(10~ 26 Wm _2 Hz _1 ) of the VC radiation in a frequency band 
Av around v, and at an angle 6 to the shower axis, from a 
charged particle shower of total energy E s inside the lunar re- 
golith can be expressed as 



F(R,6,V,E S ) = 1.89 x 10 y e 



9 „-Z 



vo[l + (v/v ) 



1.441 



sinf? 
sinf?vc 
1 Av 



(1) 



100 MHz 



Here, E = 10 20 eV, Z = («cos0 - l)/(A V c\/« 2 - 1), 
Vo = 2.5 GHz is the decoherence frequency where the wave- 
length becomes comparable to the transverse dimension of the 
shower, Rg = 1 .738 x 10 6 m the lunar radius, R the distance 
from the source point in the regolith to the detector, and Aye = 
Q.Q3Q2[vqL(Eq)]/[vL(E s )] radians the angular spread around 
0vc- The shower length L(E S ) = 12.7 + 2/3 log (E s /Eq) in 
units of radiation length 11711 . The detection threshold in eV 
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in terms of an effective antenna collection area A e g , a radio 
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FIG. 1: (color online). Monte Carlo simulation results for the frac- 
tion of neutrino events (energy E v ) above a given depth in the Moon, 
as detected by a tripole antenna (frequency v) in a Moon-orbiting 
satellite at 100 km altitude. 



system noise temperature r no i se , and a minimum detectable 
radio signal of TV 2 times the background noise. 



Based on properties of lunar soil sample returns B25I1 . we 
model the regolith down to 100 m depth as a homogeneous 
dielectric medium with a density of p = 1 .7 x 10 3 kg/m 3 and 
a radio attenuation length of / = A/(27Tntan5) m. To allow 
for uncertainties in the loss tangent, tan 8, due to metallic 
contaminants, we choose the conservative value tan 8 = 0.01 
which accounts for the available data. Within uncertainties, 
this number is also consistent with predictions made for the 
lunar bedrock layer J25I1 . Since the attenuation in both materi- 
als can be consistently described by a single value, the transi- 
tion depth becomes unimportant, and our simplified assump- 
tion of a deep, homogeneous regolith is justified for our pur- 
poses. To estimate the impact of choosing a different depth, 
below which no neutrino events can be detected, we show in 
Fig. Q]the accumulated fraction of detected events vs. depth, 
obtained from simulations as detailed below. The effective 
depth over which showers are detected is a function of both 
neutrino energy and observation frequency. The full 100 m 
depth only contributes for the lowest frequencies and high- 
est neutrino energies. Thus, the curve for v = 100 MHz, 
E v = 10 23 eV illustrates the maximum uncertainty. 

For the calculation of the emission geometry we use 0vc = 
55°, corresponding to a constant regolith dielectric permit- 
tivity e = 3 for the full radio frequency interval considered 
1 2511 . The total internal reflection angle &rm at the Moon- 
vacuum interface is complementary to the VC angle. Thus, 
for a narrow VC cone only emissions from neutrinos which 
are upward-going with respect to the local surface will escape 
the Moon. For lower frequencies, when the cone is wider, 
total internal reflection is not a significant problem since the 
interface transmission rises rapidly to its maximum value just 
a few degrees off &tir. Rays covering the A V c cone are propa- 
gated using geometrical optics, taking into account the effects 
of attenuation, refraction in the (locally flat) Moon-vacuum 
interface and internal reflection. 

In order to estimate the optimum sensitivity in some gen- 
erality, we use two different approaches to define the radio 



sensor equipment. In the first case, an isotropic antenna 
with dipole characteristics is assumed, representing a low- 
gain measurement of the complete electric field vector. This 
can be realized, for instance, using a tripole antenna 12611 . The 
antenna length is taken to be A /2 at the highest frequency. It 
will hence be electrically short over the full bandwidth, which 
ensures nearly single-mode operation. For an isotropic mea- 
surement, the noise temperature r no i se is dominated at low fre- 
quencies by the galactic background for which we use the sim- 
ple model r gal = 1.5 x 10 6 (10MHz/v) 22 K 0]. 

For higher frequencies, the predominant noise contribution 
is from the radio receiver system and the satellite, assumed 
at a nominal temperature r sys = 300 K. In the second case, a 
perfectly beam-filling antenna array is assumed, i.e. the beam 
solid angle Q. equals the solid angle of the Moon. The corre- 
sponding effective antenna area is then given by A e ff D. = X 
and the physical size of the antenna array required to achieve 
this depends on both frequency and altitude. The assumed di- 
rectivity towards the lunar surface for this case allows us to 
set r noise = r sys = 300 K. 

Unlike terrestrial measurements of RF transients, a lunar 
satellite experiment suffers no atmospherics 112811 . Likewise, 
anthropogenic noise, known to interfere badly on Earth 1 1311 . 
is favorably low in the lunar environment, in particular on 
the far side of the Moon [29]. For Gaussian thermal volt- 
age fluctuations, the single channel rate of spurious triggers, 
exceeding a level Vb, is given by T\ = 2Av x P(| V| > Vq) = 
2Av x erfc(N a /V2). At 100 MHz bandwidth, TV 2 = 25 gives 
Tj = 115 s . Requiring n-fold coincidence of independent 
measurements in a time window T, this rate can be reduced to 
r„ = T n-1 r". In practice, multiple antennas, frequency bands 
and/or polarizations are used to define the coincidence chan- 
nels. To avoid technical details in this generic study, we sim- 
ulate the detection system using an effective description with 
full rejection of thermal events and the sensitivity of a single 
channel. A threshold 7Y 2 = 25 is used, corresponding to a re- 
alistic number for each channel in coincidence (cf. |9[]). 

For an isotropic neutrino flux 4> v (E) the number of detected 
events in the energy interval [E\ ,£2]. during time At, is 



N v =At [ 2 <S> v (E)a(E)dE 
Je, 



(3) 



where the aperture function <x(E) describes the total experi- 
mental efficiency as a function of energy. Including the neu- 
trino interaction, radio wave propagation effects and the re- 
ceiver noise, it represents an equivalent detector area and ef- 
fective solid angle for which 100% of the incident neutrino 
flux is detected. This complicated aperture is best studied 
by MC simulations. For zero detected events (signal and 
background), a model independent limit 112811 on a sufficiently 
smooth neutrino flux function can be deduced by inverting (0 
over an energy interval d£ ~ E to obtain 



E<$> V {E) < 
with the Poisson factor s, 
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FIG. 2: (color online), ^--weighted flux of UHEC"v. Solid (color) 
curves show the projected detection limits from Eq. ©, based on one 
year of satellite measurements with (a) a single tripole antenna and 
(b) a beam-filling antenna for frequencies of 100 MHz (lower set of 
curves) and 1000 MHz (upper set of curves). Within each set, the 
curves from top to bottom are for satellite altitudes H of 100, 250 
and 1000 km, respectively. Dashed lines show predicted fluxes from 
the GZK process Q] (consistent with the Waxman-Bahcall bound 
5 x 10~ 8 lU), Z-bursts and Topological Defects (TD) dJ]. Thin 
solid lines show current flux limits from ANITA-lite OJ, RICE Q], 
GLUE (H] and FORTE (H]. Dotted lines show predicted sensitivi- 
ties for ANITA H, LOFAR d and LORD (H. 



To map out the most favorable conditions for UHECv de- 
tection we perform MC simulations of the aperture for several 
values of the center frequency v and altitude H for a satellite 
in stationary, circular orbit. The frequency range is limited 
to 100-1000 MHz, since the narrower VC cone strongly dis- 
favors observations closer to the decoherence frequency Vo- 
Around each center frequency, a symmetric bandwidth of ±50 
MHz is assumed. The altitudes range from 100 to 1000 km. 
For the beam-filling case, with A e s depending on altitude, also 
the distance to Earth is considered for comparison. At en- 
ergies two orders of magnitude above threshold, where the 
detection starts to become fully efficient, the apertures reach 



10 3 — 10 6 km 2 sr for moderate altitudes. 

Converting the resultant apertures to limits on the neutrino 
flux according to Eq. (01, one limit curve is obtained for each 
parameter pair (v,H). Fig. [2] shows the results for an effec- 
tive observation time of one year, the minimum duration we 
consider for a lunar satellite mission. The general trend is 
that lower frequency observations yield more stringent flux 
limits as a result of the increased angular spread Aye of the 
Askaryan radiation and the longer radio attenuation length in 
the regolith for lower frequencies. The optimum altitude is 
less clear. From geometry, the accessible aperture increases 
with altitude as the visible surface area of the Moon increases, 
giving stricter limits. However, the threshold energy also in- 
creases, which results in the successive shift towards higher 
energies of the limit curves. As shown in Fig. |2j our flux 
limits would improve substantially over the existing ones. 
They would also be competitive, and in some energy regions 
even better than, the estimated limits that may be obtained by 
ANITA |@] and LOFAR H- Predictions with both higher 
sensitivity and lower threshold can also be obtained for a lu- 
nar satellite experiment, as indicated by the LORD curve 11811 . 
assuming that more elaborate antenna setups can be realized. 

Since the threshold energy and aperture cannot be varied in- 
dependently, optimum satellite parameters can only be judged 
with respect to a specific neutrino flux model. Based on the 
fluxes from models for GZK, Z-bursts and topological defects, 
Table UJ presents integral event rates calculated from Eq. (0). 
For the Z-burst rates, we conservatively rescale the original 
flux J30ll by a constant factor to conform with the integral 
bound £ 2 4>y < 1.6 x 10~ 6 GeV cirr 2 sr" 1 s" 1 from ANITA- 
lite ISO. The rates thus obtained confirm that low frequency 
observations are most efficient. It is evident that GZK neutri- 
nos cannot be detected in the single tripole case, whereas neu- 
trinos from the other sources can. Here, the lower threshold 
energy associated with lower altitudes is strongly favorable 
for the detection of TD neutrinos, while detection of Z-burst 
neutrinos gains more from the increased aperture at higher al- 
titudes. 

In the beam-filling case the antenna gain pattern is adapted 
to the lunar disc, resulting in higher event rates, but they are 
even more affected by the choice of observation frequency. 
The suppressed galactic noise pushes the optimum frequency 
to its lowest possible value, adding to the effect of an increas- 
ing Aye- A low frequency beam-filling setup would be able to 
detect also the GZK neutrinos. For the TD and Z-burst mod- 
els, the event rates depend weakly on the altitude. 

These relatively low rates can be increased through sev- 
eral possible improvements. For CC V e events, an additional 
shower results from the emerging electron, which carries on 
average 80% of the total neutrino energy. However, for elec- 
trons of high energy (>8x 10 14 eV 113 111 ) the cross section 
to interact with the material is reduced due to the LPM ef- 
fect B2I] . The result is an elongated shower, with a reduced 
Ave 12011 . and hence smaller detection probability. Still, the 
electron-initiated showers, from the expected one third of v e 
in the primary neutrino flux, may be further investigated, as 



TABLE I: Event rates for satellite observations at frequencies v and 
altitudes H based on neutrino fluxes from the GZK process [J], a 
model for topological defects (TD, with Mx = 2 x 10 13 GeV) and 
Z-bursts (m v = 0.33 eV, consistent with cosmological data 1301 . 
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rescaled by a factor 1.5 x 10 to satisfy current flux limits). 



Satellite parameters 


Integral event rates (yr ' ) 


v (MHz) H (km) 
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TD 


Z-burst 


(a) Tripole case 








100 100 


4.9 x 10~ 2 


2.2 x 10 2 


3.3 x 10 3 


100 1000 


< 10~ 5 


0.0 


4.1 x 10 3 


1000 100 


6.8 x 10- 4 


19 


1.5 xlO 2 


1000 1000 


< 10~ 5 


1.5 x 10~' 
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(b) Beam-filling case 








100 100 


21 


6.5 xlO 3 


1.2xl0 4 


100 1000 


6.0 


1.5 xlO 4 


7.6xl0 4 


100 Earth 


1.5 


1.3 xlO 4 


1.8xl0 5 


1000 100 


2.3 x 10~ 2 


39 


1.9 xlO 2 


1000 1000 


1.2 x 10~ 3 


21 


6.2 xlO 2 


1000 Earth 


1.5 x 10~ 4 


7.4 


9.9 xlO 2 



may secondary showers from CC "v T interactions. 

In a more detailed modeling of the lunar composition, also 
the positive contributions of the bedrock layer should be con- 
sidered. The higher density means more neutrino interactions 
per unit volume and showers that develop faster, which gives 
larger Aye and increased detection probability. The effects of 
surface roughness scattering and possibly diffuse scattering at 
the rock-regolith interface should then also be considered. 

To increase the experimental sensitivity to GZK neutrinos, 
the threshold energy [Eq. ©J must be lowered by increasing 
the effective antenna area, or by decreasing N a . This means 
using more antennas in coincidence and, for the beam-filling 
case, each with a narrower beam. The sensitivity is limited by 
what is technically feasible to use on-board a lunar satellite. 

In conclusion, we have used MC simulations to demon- 
strate the prospects for a Moon-orbiting satellite to detect, via 
the Askaryan effect, ultrahigh energy cosmic neutrinos inter- 
acting in the Moon. This method can provide competitive, 
and in some respects better, conditions compared to Earth- 
based experiments. For a given lunar mission, the incremen- 
tal cost for this experiment will not be excessive. For the 
two specific model cases we consider, neutrino energies above 
10 20 eV can be covered. The flux of GZK neutrinos, which 
originate from well known Standard Model processes, can be 
discovered with suitable satellite and antenna parameters. The 
UHEC"v from proposed exotic sources can be observed even 
at a much lower flux than predicted, resulting in strict limits 
to be set or a revolutionary discovery. 
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